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Abstract. We study theoretical X-ray light curves and polarization properties of accretion-powered millisecond pulsars. We 
assume that the radiation is produced in two antipodal spots at the neutron star surface which are associated with the magnetic 
poles. We compute the angle-dependent intensity and polarization produced in an electron-scattering dominated plane-parallel 
accretion shock in the frame of the shock. The observed flux, polarization degree and polarization angle are calculated account- 
ing for special and general relativistic effects. The calculations also extended to the case of nuclear-powered millisecond pulsars 
- X-ray bursts. In this case, we consider one spot and the radiation is assumed to be produced in the atmosphere of the infinite 
Thomson optical depth. The light curves and polarization profiles show a large diversity depending on the model parameters. 
Presented results can be used as a first step to understand the observed pulse profiles of accretion- and nuclear-powered mil- 
lisecond pulsars. Future observations of the X-ray polarization will provide a valuable tool to test the geometry of the emission 
region and its physical characteristics. 
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1. Introduction 

Since the discovery of pulsars in the 1960s, determination of 
the equation of state of neutron stars based on measuring their 
mass-radius relation was one of the main research goals. Radio 
, pulsars have been used extensively for the mass determination 
(Thorsett & Chakrabartv l 1 999l l. but hardly can help to get any 
, constraints on the neutron star radius since the radio emission 
' is probably produced far from the surface. Observations of the 
isolated neutron stars emitting thermal radiation from their sur- 
face on the other hand is much more suited for the radius de- 
termination using the Stefan-Boltzmann law. Some of the radii 
determined this way from the X-ray emission (e.g. Drake et 
al. .2002 Thoma et al.lQOA-} are much smaller than any of the 
equation of state allow, pointing either to the existence of the 
more compact strange stars or (more likely) to a surface tem- 
perature gradient. 

Statistical studies of the pulse profiles from standard X-ray 
pulsars also have been used to put constraints on the geometry 
of the emission region as well as the star compactness (Bulik 
et al. l2003> . The problem here is that the accretion rate is close 
to the (local) Eddington, the magnetic field strongly affects the 
radiative transport resulting in rather complicated pulse profiles 
at low energies, and therefore it is difficult to predict from the 
first principles the angular distribution of the radiation escaping 
the emission region. 

Send offprint requests to: J. Poutanen 



Quasi-periodic oscillations (QPO) observed in a number 
of low-mass X-ray binaries (van der Klis 2000), if produced 
in the inner part of the accretion disk, could be used to put 
weak constraints on the neutron star mass as well as its ra- 
dius (see e.g. Miller et al. '1998V However, this method is in- 
direct. Discovery of millisecond coherent pulsations during X- 
ray bursts in 13 low-mass X-ray binaries (so called nuclear- 
powered millisecond pulsars, see reviews by Sti'ohmaver'2 00H 
Strohmayer & Bildsten 2003 1 and in the persistent emission 
of five sources (accretion-powered millisecond pulsars, see re- 
view by Wijnands'2007) opens a completely new range of pos- 
sibilities. The emission in these cases is produced at the surface 
of a rapidly spinning neutron star and the observed pulse pro- 
files and oscillation amplitudes are strongly affected by spe- 
cial as well as general relativistic effects (Braje et al. .2000.i 
Weinberg et al. 120011 Poutanen & GierUiiski l2003l hereafter 
PG03). Weak magnetic field in these sources does not affect 
the radiative transport and thus a much more reliable prediction 
for the radiation pattern from the surface can be obtained. This 
tightens the possible range of other parameters. However, simi- 
lar light curves still can be produced with a very different set of 
parameters. For example, exchanging inclination with the mag- 
netic inclination (angle between magnetic and rotational axis) 
has almost no effect on the light curve. 

A possible way to distinguish between the models is to ob- 
serve the behaviour of the polarization degree and polarization 
angle with phase. Polarization has proven to be a valuable tool 
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in determining the geometry of the emission region in radio 
pulsars (e.g. Blaskiewicz et al. 1991). Polarimetric observa- 
tions in the X-ray range could also be possible in the near fu- 
ture (see Costa et al. 120011 Bellazzini et al. .2003. Marshall et 
al.p003 1. In this paper, we perform theoretical calculations of 
the light curves and polarization behaviour in accretion- and 
nuclear-powered millisecond pulsars accounting for relativistic 
effects. The presented models could be used to determine the 
physical parameters of the X-ray millisecond pulsars. 



2. Method 

We assume that the emission originates in one or two antipo- 
dal (in the case of accreting ms pulsars) spots at the neu- 
tron star surface. The hard X-rays are produced by thermal 
Comptonization in a plane-parallel slab (accretion shock) of 
optical depth of order unity (PG03). The emission during the 
X-ray bursts is assumed to be produced in a semi-infinite elec- 
tron scattering dominated atmosphere. The radiation pattern is 
assumed to have azimuthal symmetry in the comoving frame 
of a spot. We compute the Stokes parameters in this frame 
are transform them to the observer frame. First, we make the 
Lorentz transformation to the non-rotating frame accounting 
for Doppler boosting and relativistic aberration and then follow 
photon trajectories to the observer at infinity in Schwarzschild 
space-time. Deviations from the Schwarzschild metric due to 
the stellar rotation have a small effect (Braje et al. 2000 1 and 
are neglected here. For pulsar rotational frequencies of v > 400 
Hz, we also need to account for time delays which in the ex- 
treme cases can reach about 5-10 per cent of the pulsar period. 

Lorentz transformation and gravitational light bending do 
not change the polarization degree and the observed value can 
be found if one knows the polar angle at which a photon is 
emitted in the spot comoving frame. For a slowly rotating star, 
the polarization vector lies in the plane formed by the spot nor- 
mal (spot radius-vector) and the line of sight as in the rotating 
vector model of Radhakrishnan & Cooke (,1969) . This vector 
rotates if a star is rapidly spinning (Ferguson ll973llT976 ). 

In the following calculations we assume two identical an- 
tipodal point-like spots, of which both, or only one is visible 
depending on the parameters. We follow the method described 
in detail in PG03 for the flux computations and extend it to the 
description of the polarization degree and polarization angle. 

2.1. General definitions 

Let 5" be a surface element (spot) at colatitude 9 (see Fig. Q] 
for the geometry). (Quantities with primes correspond to the 
frame comoving with the spot.) The second spot is situated 
at the opposite side of the star. Let / be the inclination of the 
spin axis and the unit vector in the direction to the observer 
k - (sin /, 0, cos /) lies in the (x, z) plane. The coordinates of 
the spots vary periodically and the unit vector of the spot nor- 
mal is n = ±(sin 6 cos 0, sin 9 sin (p, cos 6), where the upper sign 
corresponds to the spot closer to the observer (primary) and the 
lower sign to the antipodal spot (secondary). Thus the angle 




Fig. 1. Geometry of the problem. 



between the direction to the spot and the line of sight is given 
by 



cos ij/ = k ■ n- +(cos / cos + sin / sin 9 cos 0), 



(1) 



where the pulsar rotational phase (p is zero when the primary 
spot is closest to the observer. 

In Schwarzschild metric the photon orbits are planar and 
the original direction of the photon near the stellar surface 
is then (PG03): 



^0 = [sinff k + s,m{ijj - a) n\l sini/', 



(2) 



where a is the angle between and n, i.e. cos a - k^- n. 

The intensity of radiation and its polarization at a given pul- 
sar phase depend on the angle a' between an emitted photon 
and the local normal to the stellar surface in the spot comoving 
frame 

cos a' - 6 cos a, (3) 
where 5 is the Doppler factor. It can be expressed as 



5= l/[r(l-yScos^)], 



(4) 



where y = (1 - ff) '^^ and fi is the spot velocity in units of 
speed of light as measured in the non-rotating frame: 



2nR 



: sin 6. 



(5) 



Here the observed pulsar frequency v is corrected for the red- 
shift - rg/R, where rg = IGMjc^ is the Schwarzschild ra- 
dius of the star of mass M and radius R. The angle ^ between 
the spot velocity )8/yS = (+ sin (p, + cos (p, 0) and can be ob- 
tained from Eq. (|2j: 



cos f = — ■ kn 



sin a J3 ^ 
sinip j3 



smo' 
H sm / sm ( 



sm iff 



(6) 
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2.2. Light bending and time delay 

Due to the gravitational bending of light the angle at which 
the radiation is observed is different from the actual emission 
angle (see Fig.^. The relation between angles a and tp can be 
obtained by computing an elliptical integral (Pechenick et al. 
[T983t : 



-f 

Jr 



dr 






,2 


1,2 ^2 (l 


r 



-1/2 



where the impact parameter 
R 



■ sm a. 



(7) 



(8) 



For many purposes one can use a simpler expression derived 



(9) 



by Beloborodov ll2002l hereafter B02): 
cos a X cosi/'^l - -^j ■ 



R 



2.5rg. The spot is visible when cos a > 0. In Beloborodov's 
approximation this coiTesponds to cos i/r > -rg/(R - rg). 

The difference in light travel time around the star becomes 
larger for more compact stars. A photon of impact parameter b 
is lagging one with the impact parameter by (Pechenick et al. 



cM{b) 



dr 



1 - 



-1/2 



-1 



(10) 



For rough estimations one can use an approximation At{b) » 
R{\ — cos (/')/c. If we compute the phase delays relative to pho- 
tons emitted at = 0, then the observed phase is approximately 

<^obs ~ + A0o(l - COS0), (11) 

where 

A0() - cos / cos 6 InvR/c . (12) 

For small delays, this expression can be reversed to 

cos 4> X COS 0obs + A0O sin 0obs(l - cos 0obs), (13) 

Sin0 a: Sin0obs + A0o cos 0obs(l - cos 0obs)- 



2.3. Observed intensity and flux 

The combined effect of the gravitational redshift and Doppler 
effect results in the following relation between the monochro- 



matic observed and local intensities (see e.g. Misner et al. ll973l 
Rybicki & Lightman fDT^ : 



where E/E' - 6 ^Jl - r^/R. Here I'^,{a') is the intensity com- 
puted in the frame comoving with the spot. For the bolometric 
intensity, one gets 



{6^\-rg/R^ I'ia'). 



(15) 



If the radiation spectrum can be represented by a power-law 
I'pia') oc E'^^^^^^ (as in accretion-powered ms pulsars) with the 
photon spectral index F which does not depend on the angle a', 
then 



I',,(a')^I',(a'){6^1-rg/Rf \ 



(16) 



This approximation is equivalent to the assumption of a weak 
energy dependence of the angular distribution. In the case 
of thermal Comptonization, it can be used if the maximum 
Doppler shift AS ~ InvR/c sin/ sin 6 is smaller than the 
typical (relative) energy separation between photons of sub- 
sequent scattering orders AE/E ~ AkT^-lm^c^, where 
is the electron temperature. This condition is equivalent to 
(v/600 Hz) sin/ sinfl < kTf./ 16 keV, which is always sat- 
isfied in accreting pulsars. Then, the observed photons at a 
given energy E belong mostly to a certain scattering order 
(n = 0, 1, 2, ...) which can be determined from 



which gives an accuracy of the order of a few percent for R > 



(17) 



where £0 is the mean energy (at infinity) of the seed photon 
black body distribution. 

In order to compute the observed flux, we need to know 
variations of the solid angle occupied by the spot on the ob- 
server's sky (PG03): 



S ' cos a' 



1 



dcosff 



D2 



1 - r„/R dcos tfr 



(18) 



where D is the source distance. (This transforms to 
Q a! S'cosa'/D^ in Beloborodov's approximation, since 
dcosQ'/dcost/' a 1 - rg/R.) Hence the monochromatic flux 
becomes: 

Fe = /^Q = (1 - rg/R)''^6%,(a')cosa^^^, (19) 

d cos (f/ 

where we used aberration formula (|3j- The bolometric flux is 
given by: 

dcos a S' 



F ^{l- rg/R) 6' I' (a') cos a- 

a cos iff 

For a power-law spectrum, the monochromatic flux is 

r^t , , dcOSCK^' 

Fe^(\- rglRf^^^'^ 5^^^ I'Eia) COS a 



(20) 



(21) 



d cos i]/ 

In the following calculations we have used Eqs. MQ\ and MW 
(with F = 2) for the black body and power-law spectra, 
respectively. We normalize the fluxes with Fq - 2(1 - 



(14) 2.4. Polarization degree and polarization angle 



Relativistic effects do not change the polarization degree and 
thus we take it equal to the polarization in the spot comoving 
frame P'p(a'). However, when calculating the polarization de- 
gree of the total flux from both spots, the difference in the tilt 
of the polarization ellipses (difference caused by the time delay 
and different relativistic rotations) need to be considered. 
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In order to describe polarization, we introduce the main po- 
larization basis. 



0} - cos / k 



sin I 



kXOJ 



smi 



(22) 



where d) denotes the unit vector along the stellar rotational axis. 
In the absence of relativistic rotation of the polarization plane 
(for a slowly rotating star), the polarization vector lies in the 
plane formed by the local normal n and the direction to the 
observer k. The corresponding polarization basis is 

(23) 



n - cos {jj k 

e\ - : ; , 62 



sin ip sin ip 

The polarization angle {Y'A)xo measured from the projection of 
the spin axis on the plane of the sky in the counter-clockwise 
direction is given by: 



COS^O = -^1=^2 ■ ^2 = 

sin^o = ■ ^1 = -ef ■ ei 
We thus get 



sm; cos 6/ -cos; smt/ cos( 



sin)/' 



sin0 sin^ 
sini/r 



tan;fo 



sin0 sine 



sm I cos y — cos ; sm a cos ( 



(24) 



(25) 



This expression neglects the rotation of the polarization plane 
due to relativistic motion (e.g. Ferguson Wl^^WW). To cor- 
rect for that we introduce the polarization basis related to the 
photon direction close to the neutron star surface 



n — cos a ko o 
smc 



^0 X n 

2 - —■ = ^2- 

sma 



(26) 



The equality = C2 is related to the fact that photon tra- 
jectories are planar. If the polarization vector in the spot rest 
frame lies in the meridional plane, then it will be transformed 
to (Nagirner & Poutanen[T993} 



, n + 6y cos a (B - ko) 
By — 

Vl - <5^ cos^ a 



(27) 



in the non-rotating frame. Thus, the rotation of the polarization 
vector due to relativistic motion is given by 



cos;t'c - e\ ■ e\ , sin;^c - e\ ■ e\ , 
which is reduced to (Poutanen, in preparation) 
cos / sin - sin / cos 6 cos S 



(28) 



tanxc - jScos a- 



sin a sin iff + p sin / sin ^ 



This correction angle depends linearly on the spot velocity /3 
and, therefore, is negligible for spin frequencies v < 300 Hz. 
Since the angle Xc is constant for parallel transport along the 
photon trajectory, the total polarization angle for each spot is 



X^Xo +Xc 



(30) 



When we plot polarization angles for individual spots, we ro- 
tate it by 7t/2 (in any direction) if the polarization P - Q/I is 
negative (i.e. dominant electric vector oscillations are perpen- 
dicular to the meridional plane defined by the spot normal and 
the direction of photon propagation). 



For each spot calculations provide us with the Stokes vec- 
tor (F/, Fq, 0)^ (we do not consider sources of circular polar- 
ization), where flux Fj is computed using Eqs. j20> or (12 1> . 
Fq = PFj, and the polarization degree P at angle a' is com- 
puted for a given optical depth of the slab (see below Eqs. J39> 
and J41» . We rotate this vector to the main basis 



1 
cos 2x - sin 2x 
sinlx cos2;t^ 





( Fj ] 




' Fi ' 




Fq 




Fq cos 2x 




[ J 




. Fq sin2x , 



(31) 



Denoting the corresponding quantities for the primary and sec- 
ondary spots by indicies p and s, respectively, we get for the 
total observed Stokes vector 



F^ 



F^Q cos 2xp + F'q cos 2xs 
sin2^p + F^sin2^s . 

The degree of polarization is then obviously 



(32) 



■iFl>)^- 



■ 2F'^qF'q cos(2^p - 2x,) 



r J -r r J 



and the total polarization angle is given by 



tan2^tc 



Fq sin2^p h-F^ sin2;\-s 
F^Q cos 2xp + F'q cos 2xs 



(33) 



(34) 



3. Local intensity and polarization 

3.1. Radiative transfer in the accretion shock 

We assume that the accretion shock can be represented as a 
plane-parallel slab which horizontal size is much larger than 
its vertical height so that the photons escaping from the sides 
of the slab can be ignored. This is a good approximation for 
moderate luminosities (Basko & Sunyaev 1 19761 typical for 
accretion-powered millisecond pulsars (Wijnands |2004< . We 
calculate the angle-dependent intensity and the polarization de- 
gree of radiation emerging from the emission region, assuming 
Thomson scattering and coiTesponding polarization. The cal- 
culations are carried out following the procedure described in 
Sunyaev & Titarchuk (1985 1. Unpolarized seed soft photons 
are assumed to be injected from the neutron star surface (i.e. 
bottom of the slab). The boundary conditions are 



(29) hr(T = 0,/^) = I,,At = to, -/^) = 0, < // = COS < 1, (35) 



where 6 is the angle between photon direction and the outward 
slab normal, indices /, r refer to the polarizations in the merid- 
ional plane and perpendicular to it, respectively. The Thomson 
optical depth t is measured from the bottom of the slab. The 
intensity of unscattered photons is given by: 



(36) 



The intensity of n times scattered photons can be computed 
from 



/r(T,yu) 



,3, 



K. Viironen and J. Poutanen: Light curves and polarization from millisecond pulsars 
I ' ' ' I ' ' ' I ' ' ' d 60 



40 - 




- Q, 



1 1 1 1 


1 1 1 1 1 1 1 1 1 

30 


1 1 1 1 1 1 


- 




- 




3 






2 






8 

n— 

1 1 1 1 1 1 1 1 L 









2.0 



1.5 



1.0 



0.5 - 






- a. 



10 




1 30 

■ — "10 : 



1 1 1 1 1 1 















-10 




5^ 




-20 




-""^^ 2/^ 




-30 








-40 




Tl=l/ 




-50 




1 , , , 1 , , , 1 


, , 1 , , , 



0.2 



0.8 



1.0 



0.0 0.2 



0.4 0.6 



0.8 



1.0 



Fig. 2. Intensity (normalized to unity at /u - 1) and polarization of the radiation escaping from a slab of Thomson optical depth 
To = 0.5, 1,3 for different scattering orders n. The unpolarized seed photons with intensity = 1 are injected from the bottom 
of the slab. The dashed curves show the intensity and polarization corresponding to the classical results of Chandraselchar-Sobolev 
for the semi-infinite atmosphere tq = oo. 



{(1-//W')+//' [B"(t') + C"(t')]} 



dr' 



(38) 



I';(r,^,) = j|%xp(-l-i:)[B«(r') + C«(T')]^, (37) c"(r) = [/^'(r,/) + /^'(r, -/)] d;. 

where the quantities A"(t), B"{t) and C"(t) are as follows 
A"(r) [tr\r,fi') + ir\r, -fi')] d/. 



The intensities /{'^(t, -//) can be computed in a similar way by 
changing the integration boundaries in Eqs. i37i from (0, t) 
to (r. To). These equations are similar to those in Sunyaev & 
Titarchuk ( 1985 1, but we do not assume symmetrical (relative 
to the slab central plane) distribution of seed photons. Finally, 
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the degree of polarization of the radiation leaving the slab up- 
per boundary is 



p^Q{To,^i)|I{To,^l), 

where the Stokes parameters are 

/(ro,jU) = Il(To,fl) + Ir(T(hfl), 
g(T0,//) = Il(To,fJ.)- IriTQ,fJ.). 



(39) 



(40) 



The results of calculations for different scattering orders are 
presented in Fig.|2l 

Due to the symmetry, the polarization vector can be either 
in the meridional plane (positive polarization) or perpendicular 
to it (negative polarization). The polarization degree is zero for 
the seed radiation (n = 0), negative for small scattering orders 
and for tq = 1 becomes positive at n > 3. The polarization 
degree also depends on the electron temperature. At high elec- 
tron temperatures, the electron random motions produce dif- 
ferent aberrations and in the electron rest frame photons are 
scattered at different angles. This reduces the polarization de- 
gree (Poutanen 1994; Poutanen & Svensson 1996). Thus, for 
example, for kT^ ~ 100 keV, polarization is two times smaller 
than that given in the Thomson scattering approximation. 

Variation of the optical depth also strongly affects the an- 
gular distribution of radiation and the polarization degree. An 
increase in tq leads to more beaming along the normal to the 
slab and to smaller P (Sunyaev & Titarchuk 1985). At larger 
T() ~ 3 and n » 1 (see n = 10 curve at the lower panel in 
Fig. O the polarization changes the sign for some angles re- 
minding the behavior of P for tq = 1 and n - 3. 

Since in the process of thermal Comptonization the energy 
of a photon increases with the number of scattering, there is 
one-to-one correspondence between scattering order and the 
energy. For our further calculations we choose optical depth 
To = 1 since the X-ray spectra of accreting neutron stars at 
low and intermediate luminosities are often described by ther- 
mal Comptonization in a slab of optical depth of order unity 
and electron temperature kTf. ~ 50 keV (Gilfanov et al. i ll 9981 
Barret et al. 2000 PG03). Polarization degree for larger and 
smaller tq can be estimated from Fig.|2] 

3.2. X-ray bursts 

If during the X-ray bursts, we can observe coherent pulsa- 
tions, the atmosphere is pinned down to the neutron star sur- 
face, and thus a plane-parallel approximation is still valid. In 
a classical problem of radiative transfer in a plane-parallel 
semi-infinite (tq — > oo) atmosphere with opacity domi- 



nated by Thomson scattering (Chandrasekhar & Breen ll947l 
Chandi-asekha nl960t Sobolev il949. Sobolev 1963^, the in- 
tensity of escaping radiation is approximated by the formula 
/(/i) ~ 1 + 2.06yU. The degree of polarization is at its maxi- 
mum, 1 1 .7 %, at ju = and decreases to 0% at // = 1 due to the 
symmetry. The polarization degree can be approximated as 



■11.71' 



(41) 



1 + 3.582yL< 

The dominant direction of the electric vector oscillations is per- 
pendicular to the meridional plane. 
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Fig. 3. Classes of the light curves for two antipodal spots and 
R = 2.5rg according to Beloborodov (2002). Asterisk shows 
our basic parameter set and circles correspond to other consid- 
ered cases. 



4. Results 

4.1. Pulse profiles for black body spots 

The light curves produced by two small antipodal spots emit- 
ting as a black body can be divided into four classes (B02, see 
Fig.O: I - only one spot is visible, II - the primary spot is al- 
ways visible, the secondary spot is visible sometimes. III - both 
spots are eclipsed sometimes, IV - both spot are visible all the 
time. The light curve from a single spot in Beloborodov's ap- 
proximation can be represented as follows (PG03; Poutanen & 
Beloborodov, in preparation): 

F oc 6^[rg/R + (1 - rg/R)cosiJ/] 



^1-5/3 ^1 - rg/R sin / sin (2 -H t/ cos (p), (42) 



where 

U = (1 - rg/R)smi sind, 

Q - rg/R + (l - rg/R) cos / cos 0. 



(43) 



Here we took into account only the first term of the expan- 
sion of the Doppler factor in velocity /3 and approximated 
sin a/ sin by its asymptotic value for small angles ^Jl - rg/R. 
Thus for small /3, the observed flux varies as a sinusoid. At high 
rotational frequencies, the oscillations are shifted in phase and 
a harmonic appears as a result of Doppler effect. An additional 
source of the phase shift and the signal at the harmonic is the 
light travel delays (see Eq. il3\ ). One can notice that the ex- 
change of / and does not affect the light curve (for small 
/3) since the flux depends on the combinations sin / sin and 
cos i cos (see Eqs. S)- 
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Fig. 4. Observed flux from two antipodal black body spots as a function of phase 0obs/(27r). The basic set of parameters is6 - 20°, 
/ - 60°, V = 300 Hz, R = 2.5rg and M = 1.4M0. a), b) Dependence on the inclination angle / and the magnetic angle 6 with 
other parameters from the basic set. Dependence on the neutron star compactness is shown in panels (c), while panels (d) show 
the effect of different spin frequency. Solid curves correspond to the total flux, while the dashed curves are for the primary spot, 
and dotted curves - for the secondary. 



In Fig.|^,b we plot the light curves for different couples of 
i, 6 shown in Fig.|3l The phase shift due to the light travel effect 
is computed relative to the photons emitted at <p = (so that 
the observed phase coincides with the pulsar phase </> = 0). 
For a given 0, the variability amplitude from a single spot is 
increasing with / since polar angle a covers a larger interval. 
The semi-amplitude is simply 

A = - F„i„)/(F„,ax + i^min) ^ U/Q. (44) 

The three panels on the top raw show the effect of vary- 
ing inclination only. At small inclinations, only primary spot 
is visible (class I) and deviations from the sine wave are small. 



At larger inclination, the second spot appears at some moments 
(class II), when the flux from the first spot is at minimum. The 
total variabiUty ampUtude then decreases. For a slowly rotating 
star, the increase in the second spot flux would exactly cancel 
the decrease of the flux from the first spot making a plateau 
in the total flux. This does not happen when the star rotates at 
300 Hz. Rapid rotation introduces an asymmetry, boosting the 
radiation when the spot is moving towards the observer. Since 
the peak of the Doppler factor reaches the maximum quarter of 
the period earlier than the projected area, the peak of the light 
curve shifts to an earlier phase. At very large inclinations (third 
panel a), both spots are visible all the time. The Ught curve be- 
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Fig. 5. Pulse profiles of the flux, polarization degree, and polarization angle of the scattered radiation for different scattering 
orders. The emission region is a slab of Thomson optical depth tq = 1 . The same basic set of parameters and the same designations 
for the curves are used as in Fig.|4] 



longs to the class IV and in the absence of rotation would be 
flat (B02). The variability amplitudes at the fundamental fre- 
quency and its harmonic, proportional to oc v sin 2i sin 29 and 
to oc V sin^ / sin^ 9, respectively (Poutanen & Beloborodov, in 
preparation), are small. 

The light curves for other values of inclination / and the 
magnetic angle 9 are shown in Fig.|^. One sees that exchange 
i and 9 has no effect on the light curve as was noticed above 
on theoretical grounds (see Eqs. J42> . J43l l). The variability 
amplitude for individual spots increases with 9 since it is ap- 
proximately proportional to sin 9 as can be seen from Eq. ( I44> . 
At very large / and 9 (rightmost panel Fig. |^), both spots 
are sometimes eclipsed. Two plateaus, that would appear for 
slowly rotating stars when both spots are visible, are not present 
when rotation is fast. Two pulses corresponding to the two spots 
are almost identical producing a strong signal at a harmonic 
frequency. 

Figure 15 demonstrates the dependence of the light curves 
on the neutron star compactness MjR. For large compactness, 
i.e. small radius R < 2rg, both spots are always visible (class 
IV) for most of the parameters /, 9. The variability amplitude is 
a linear function of the spot velocity. For larger radii, the second 
spot becomes eclipsed (for the set of parameters chosen), and 
the variability amplitudes grows. This can be easily understood 
from Eq. i44i . 



The dependence on the rotational frequency is shown in 
Fig-ISl. For very slow rotation speed the individual profiles are 
close to pure sinusoids and the total profile is symmetric with 
flat parts at the phases were both spots are visible (B02). The 
larger is the spin frequency, the more the individual profiles 
become skewed and the more distorted is the total profile. The 
individual pulse amplitudes also slightly increase with increas- 
ing frequency due to variations in the Doppler factor. 

4.2. Pulse profiles and polarization of the 
Comptonized emission 

The spectra of accreting millisecond pulsars can be represented 
as a sum of a black body-like emission and a Comptonized tail 
(Gilfanov et al. [T998l Gierlinski et al. l2002l PG03; Markwardt 
et al. 2002' Galloway et al. l2002t . The black body photons play 
also a role of seed photons for Comptonization. Since below a 
few keV contribution of the black body is large the expected 
polarization degree is small (unless this emission is a Wien 
peak from electron-scattering dominating atmosphere which 
can be polarized). Scattering in the hot electron slab modifies 
significantly the angular distribution of radiation and produces 
polarization signal as described in § |3] (see also Poutanen & 
Svensson[1996l. The polarization degree is a strong function 
of the scattering order for small n (see Fig. |2j and therefore 
the energy (close to the peak of the black body). After 7-10 
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Fig. 6. The amplitudes of the fundamental fli(n) and the first harmonic a2{n) (thin curves, fundamental has a larger amplitude) 
and the total rms (thick curves) as a function of the scattering order are shown in the top panels. Solid, dotted, and dashed curves 
correspond to the slab optical depth tq = 0.5, 1,3, respectively. The model phase lags at the fundamental frequency relative to 
n = 0, <p\(n) - <p\(Q) (thick curves), and those at the first harmonic, 02(n) - 02(0) (thin curves), are shown in the lower panels. Left 
panels (a) correspond to one spot and right panels (b) to two antipodal spots. The parameters are: / = 60°, = 20°, R - 2.5rg, 
M - \ .4M0. The pulsar frequency was chosen v = 401 Hz which corresponds to that of SAX J1808. 4-3658. The observed phase 
lags from that source are shown in lower panel (a) by triangles. We transformed the photon energy to the scattering order using 
Eq. Mil with kT^ - 50 keV and £0 - 2.8 keV. Diamonds represent the phase lags observed in the 2.3 ms (435 Hz) accreting 
pulsar XTE J1751-305 (Gierlifiski & Poutanen, in preparation). 



scatterings polarization saturates. For most of the simulations 
below we use n - 7. This would correspond to about 10 keV for 
a typical electron temperature of 50 keV and the seed photon 
temperature of 0.6 keV. At higher T^, the corresponding energy 
is higher, while polarization at 10 keV is then smaller Increase 
of the electron temperature generally have a depolarizing ef- 
fect (Poutanen il994i) . We choose the slab optical depth tq - 1 
which is consistent with the spectra of SAX J1808.4-3658 
(PG03). The predicted polarization degree anti-correlates with 
the value of tq (a smaller tq would increase it and vice versa). 

Since the polarized flux coiTesponds to a power-law like 
Comptonized spectrum, we use Eq. i21\ for the calculations. 
The polarization degree for individual spots is determined by 
the polar angle a' in the spot comoving frame and we plot 
its absolute value. The total polarization degree is given by 
Eq. (l33t . The polarization angle is given by Eqs. (I25t-(l30l. If 
the polarization degree P(a') is negative (in our convention this 
corresponds to the electric vector oscillations perpendicular to 



the meridional plane formed by a photon momentum and a lo- 
cal normal to the slab), we rotate the PA for individual spots by 
90°. The total polarization angle is defined by Eq. i34\ . 

First, in Fig.|5]we show the dependence of the flux and po- 
larization profiles on the scattering order Radiation scattered 
once is beamed along the normal to the slab similar to a black 
body radiation (/(yu) ~ const). The flux profiles are thus almost 
identical to those for the black body spots (compare left panel 
in Fig.|5]to the 2nd panel of Fig.^) with the maximum close 
to the zero phase. The polarization degree is always negative at 
first two scattering orders as can be seen from Fig.|5] The PA is 
then varies around 90°. The polarization degree is larger for the 
secondary spot because we see it at grazing angles. However, 
the observed flux is larger for the primary spot and hence the 
total degree of polarization is determined by the primary. For 
each spot, the polarization is larger when the flux is lower be- 
cause polarization is a monotonic function of jj - cos a' which 
also controls the flux. 
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Fig. 7. Same as Fig.|5] but for scattering order n = 7 and varying inclination i and magnetic angle 0. 
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Fig. 8. Same as Fig.0 for another set of parameters and /. 



With increasing n the escaping flux becomes less dependent 
on the polar angle resulting in the reduction of the variability 
for a primary spot. The intensity at large angles from the nor- 
mal increases causing a relative increase of the flux from the 
secondary spot (which is visible at large angles). The total flux 
is now peaked at phase 0.5 when the secondary has the max- 



imum. At the third scattering, the polarization degree defined 
by Eq. j39t changes sign for smaller polar angles, being still 
negative for larger ones. The PA for the primary spot (which is 
observed closer to face-on) now varies around 0°, while the sec- 
ondary spots, observed at grazing angles, demonstrates a jump 
in the PA when polarization changes the sign. For scattering 
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Fig. 9. Same as Fig.0 but for varying neutron star compactnesses R/rg. 



order n - 1 the polarization degree defined by Eq. ( I39> is al- 
ways positive, and thus, the PA for both spots varies around 
zero degrees. The polarization degree is now not a monotonic 
function of //, but reaches a maximum at yU ~ 0.2 (a' ~ 80°). 
One sees that for the secondary spot P does not have the maxi- 
mum when the flux is minimal since we see the spot at grazing 
angles (see Fig.|3. The hardly visible deviations in PA of the 
individual spots are due to a combined effect of the time delay 
and a different relativistic correction angle given by Eq. j29l l. 

Let us represent the light curve for a given scattering order 
n as a Fourier series 

F„((p) = ^ |l + X "■'■^"^ cos[27r./(0 - 0/n))]| , (45) 

where F„ is the mean flux. For a pure sine wave, the rms is 
ai/ V2. The dependences of the rms and the variability ampli- 
tudes fli, 122 on the scattering order n are shown in Fig.|6lfor 
three optical depths tq. The amplitude of the fundamental is 
generally a decreasing function of n for one spot (Fig.|6^). If 
two antipodal spots are present (Fig. (SJi), the rms has a min- 
imum at n = 1-2 when tq - 0.5 - 1. At n = 0, the black 
body is peaked as /(ju) = Q-'^o/f along the normal to the slab 
producing large variability. Radiation scattered once has now a 
radiation pattern which actually is close to a pure black body 
not covered by a hot slab (see Fig. |2] middle panel). This is 
known to reduce the variability (B02, see Fig. |5^). At larger 
n, radiation becomes beamed along the slab surface, secondary 
spot becomes stronger and variability increases (Fig.|5j;). 

The phase lags ^/n) - (pj(0) relative to unscattered radi- 
ation, i.e. n = 0, at the fundamental frequency are shown in 
lower panels of Fig.|6l The negative lags mean that the harder 



photons are leading the softer ones. With increasing scatter- 
ing order the pulse profile shifts to an earlier phase. For one 
spot (left panel), the phase lag results from a combination of 
Doppler boosting and a change in the emission pattern. The 
phase lag is then a smooth function of n. It saturates after 5-10 
scattering, when the distribution of photons over the slab does 
not change anymore. The lags are larger for smaller tq when 
radiation is stronger beamed along the slab (fan-like emission 
pattern, see top panel in Fig.O. 

For two spots (right panel), there is a jump by almost half- 
a-period in the lag dependence at small n and tq. This is again 
a result of dramatic change in the radiation pattern at small n. 
At « = 0, the light curve is dominated by the primary spot and 
the peak is close to = 0, while at larger n, the peak is reached 
at phases where the secondary spot is visible. For tq = 3, the 
radiation pattern changes gradually with n and the phase lags 
also show a gradual evolution. 

Figures0|Slrepresent the light curves and polarization pro- 
files for n = 7 and different couples of / and shown in Fig.|3l 
We can see the dependence on the inclination at the first three 
panels of both figures. When / and are small, we see only 
one spot (class I) and the variability amplitude oc sin/ sm0 
is also very small. Actually, it is much smaller than that for 
the black body spot since the flux escaping from the slab af- 
ter a number of scatterings is more isotropic. The polariza- 
tion degree is small since we look at the spot almost along 
the normal. If / - 0, then according to Eq. ( I25> . the PA can 
be presented as tan;^' x -l/[cos/ tan(0/2)]. For small /, this 
reduces to a linear relation x ~ 0/2 - 7r/2 clearly visible in 
Fig.Qi; while for larger / the x-ip relation is more complicated 
(see Fig. |SJ),d). For small /?, the PA crosses the zero at phase 
4> * -2yS/( ^1 -rg/Rsin i). 
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Fig. 10. Same as Fig.0 but for varying spin frequencies - 



With increasing inclination, the second spot becomes visi- 
ble and the oscillation amplitude grows. In class IV, it is actu- 
ally larger than that for black body spots (compare Fig. Q; to 
the 3rd panel in Fig. |^) since the variations of the flux from 
the two spots do not exactly cancel each other The polariza- 
tion also grows with / since we see the spots are larger angles. 
In Fig. 05 we also see that polarization from the primary spot 
anti-correlates with corresponding flux, while for the secondary 
spot they are correlated. This is related to the fact that the pri- 
mary spot is observed at small polar angles jj ~ 1 where flux 
and polarization (see Fig.|5J are anti-correlated, while the sec- 
ondary spot is observed at larger angles where the relation re- 
verses. At larger inclinations (see Figs. Q;, |SJ),c) the relation 
between the polarization (from the secondary) and the flux is 
not monotonic since we observe the spot at intermediate an- 
gles. The PA varies around zero and is close to a sinusoid 
X ~ -(tan 6*/ sin/)sin(/> at large inclination and small 6 (see 
Figs. 03,c,|Sl;). When both / and 9 are large (e.g. Fig.lSJl), the 
flux and polarization degree vary mostly at the harmonic of the 
spin frequency since the visibility conditions for both spots are 
similar. The PA, however, varies with the spin frequency. 

As we pointed out in § 14.11 the flux profile is (almost) in- 
variant to exchanging / and 6. The degree of polarization is also 
invariant. The PA, on the other hand, show completely differ- 
ent behavior: when / > the variations are small, but in the 
opposite case PA make a full turn of 360° (compare Fig. 05 to 
Fig. Ell and Fig. |7l; to Fig. ITjl). At / < 0, the PA is almost a 
linear function of the phase, with the derivatives 



dx 
d4> 



sin dx 
^0 sm{0 - i) ' dip 



sin 9 
sm(0 + i) ' 



(46) 



The compactness of the neutron star influences the gravita- 
tional bending and changes the angle the spots are observed at. 
For small radius, bending is strong and both spots can be visible 
all the time (class IV light curve, see Fig.|5^), while at larger 
radius the light curve class can change to II or even I (Fig.|5j3- 
d). The polarization degree slightly increases with R (effect is 
clearly visible for the primary spot), since for the given angle 
1/', cos a X rg/R + {l -rg/R) cos is a monotonically decreasing 
function of R. The effect of the radius on the PA is tiny since 
only (small) correction angle i29\ depends on it. 

The influence of the neutron star spin is shown in Fig. 1101 
Increasing v modifies mostly the flux profile, since it is pro- 
portional to the Doppler factor in power F H- 3. The peak shifts 
towards the phase where the Doppler factor reaches the maxi- 
mum (a quarter of the period before the peak of the projected 
area). The pulse profiles thus becomes more skewed. The polar- 
ization degree depends on cos a' - 6 cos a, which varies only 
slightly with the spin. The PA is affected through the correction 
angle which is a linear function of v (see Eq. j29H and is small. 
The time delay slightly shifts all the profiles, especially those 
for the secondary spot. 

4.3. Phase lags in accreting millisecond pulsars 

The detailed fitting of the light curves predicted by our mod- 
els to the data is outside the scope of this paper (see e.g. PG03 
for applications to SAX J1808.4-3658). We just compare here 
the general evolution of the pulse profile with energy which 
is quantified by the phase lags. We show the phase lags ob- 
served from the first accreting ms pulsar SAX J1808. 4-3658 
(Cui et al. T998"; Gierlinski et al. "^OOTl in the lower panel of 
Fig. lit. One sees that the observed behavior is consistent with 
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Fig. 11. Light curves, polarization degree and polarization angle expected from an X-ray burst. Semi-infinite electron scattering 
atmosphere is considered. The parameters are M = I.4M0, / = 60°, R - 2.5rj, v - 400 Hz. 



the model predictions for a single spot and optical depth of 
the slab of tq = 1. Such an optical depth is in turn consis- 
tent with the observed X-ray spectra from the source (PG03). 
SAX J1808.4-3658 shows almost sinusoidal variations which 
cannot be described by two-spot model. As is shown in PG03, 
most probably, the second antipodal spot is blocked from the 
observer by the accretion disk which inner radius has to be 
smaller than a few neutron star radii. 

The second accreting ms pulsar XTE J 175 1-305 has 
pulse profiles which are almost identical to those of 
SAX J1808. 4-3658 with somewhat smaller rms amplitude 
(Markwardt et al. '2002 ). Again, only one spot seems to be vis- 
ible. The spectral analysis (Gierliriski & Poutanen, in prepara- 
tion) shows that the optical depth of the Comptonizing medium 
is To ~ 2. According to the model, the phase lags then are 
smaller than those for to ~ 1 . Observations confirm (diamonds 
at the lower panel Fig. (6^) that this indeed is the case. One 
should note here that XTE J175 1-305 has a slightly different 
pulsar period (2.3 ms vs 2.5 ms in SAX J1808.4-3658), there- 
fore the model predictions should be changed accordingly. 

We emphasize that the model phase lags are computed rel- 
ative to the zeroth scattering order (black body), while the ob- 
served lags are computed relative to the ~ 3 keV energy band 
where photons scattered in the hot slab can provide significant 
contribution. It is encouraging that this simplified model pro- 
vides a good description of the data. At this this stage, however, 
a detailed comparison seems unjustified. 



4.4. X-ray bursts 

The X-ray bursts often show coherent oscillations (see 
Strohmayer & Bildsten i2003> . The energy dissipation takes 
place deep in the electron scattering atmosphere, and thus we 
can assume that Thomson optical depth is infinite. The escap- 
ing intensity in that case is I{fi) cc I + l/u and polarization 
is 11.7% at maximum when the spot is viewed edge-on (see 
dashed curves in Fig. |2}. The predicted pulse profiles and be- 
havior of polarization degree and angle are shown in Fig. 1111 
Since the radiation pattern is peaked sharper along the normal 
than that of the black body, the variability amplitude is larger 
(compare Fig.ll Ih to 2nd panels in Fig.|4t,b). The polarization 
is increasing with the spot colatitude reaching the maximum 
of ~ 12% close to the eclipses. The PA varies around 90° as 
the electric vector is predominantly perpendicular to the merid- 
ional plane and the amplitude grows with 6. 



5. Summary 

In this work we have developed a formalism to compute wave- 
forms and polarization profiles from accreting millisecond pul- 
sars and coherent oscillations in X-ray bursts. The radiation 
from emission region was assumed to be produced as a results 
of Compton scattering in the accretion shock of Thomson opti- 
cal thickness ~ 1 (in the case of accretion-powered sources) 
or optically thick electron atmosphere (in the case of X-ray 
bursts). The observed flux and polarization are affected by spe- 
cial and general relativistic effects, such as Doppler boosting, 
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relativistic rotation of the polarization plane, light bending, and 
gravitational redshift. The developed technique is accurate and 
extremely efficient from the computational point of view. 

Future polarimetric observations in the X-ray domain will 
serve as a powerful tool in determining the geometry of the 
emission region in rapidly rotating neutron stars showing co- 
herent millisecond oscillations. Our simulations are the first 
step in theoretical modelling of these sources. 
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